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Abstract 

^\ • This paper is devoted to the problem of error detection with quantum codes. In the 

first part we examine possible problem settings for quantum error detection. Our goal is to 
derive a functional that describes the probability of undetected error under natural physical 
assumptions concerning transmission with error detection with quantum codes. We discuss 
possible transmission protocols with stabilizer and unrestricted quantum codes. The set 
of results proved in part I shows that in all the cases considered the average probability 
of undetected error for a given code is essentially given by one and the same function of 
its weight enumerators. This enables us to give a consistent definition of the undetected 
error event. In part II we derive bounds on the probability of undetected error for quantum 

V^ ', codes. 

In the final section of the paper we examine polynomial invariants of quantum codes and 
show that Rains's "unitary weight enumerators" jl6| are known for classical codes under 
the name of binomial moments of the distance distribution. As in the classical situation, 
they provide an alternative expression for the probability of undetected error. 
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£3 . 1 Introduction 

cr 

The possibility of correcting decoherence errors in entangled states was discovered by Shor in 
p9| and Steane |20|]. Since then the theory of quantum codes has been a topic of intense study. 
Error processes in the depolarizing channel are characterized in [jll]] where it was shown that 
one can restrict attention to error operators given by Kronecker products of Pauli matrices. 
This opened the possibility of finding parallels between the theory of classical error correcting 
codes and their quantum counterparts. 

This theory has developed simultaneously in different, often not very related directions. 
Mathematically one of the most interesting works on quantum codes is the paper by Calder- 
bank et al. H, where the authors introduce a class of so-called stabilizer quantum codes 
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(independently found in [[?]]), and discovered a beautiful geometric and group-theoretic con- 
nection that shows that such codes can be studied using the well-developed theory of classical 
codes over GF{4). Stabilizer codes share some of the properties of classical linear codes. 

One of the most important concepts in classical coding theory is the notion of decoding, 
decoding error, and the probability of this event for a given (random) error process in the 
channel. In the quantum case a similar theory still awaits construction; even defining decoding 
does not seem an easy problem. In the classical coding theory the simplest known decoding 
algorithm is testing the received vector for containment in the code; if the test fails, the decoder 
detects an error. Clearly, the error will not be detected if the sum of the transmitted vector 
and the error is itself a code vector. 

The focus of our paper is error detection by quantum stabilizer and nonstabilizer codes. To 
give an analogy, let us recall the definition of error detection for classical codes. Let D be a code 
of length n over an alphabet of q elements and Bi, < i < n, its Hamming distance distribution 
given by Bi = jjyrJ^c'eD \{ c e D I dist(c,c') = i}\. Suppose that D is used for transmission 
over a q-aiy memoryless symmetric channel in which each symbol is unchanged with probability 
1 — p and replaced by another symbol with one and the same probability -^r- Let c G D be a 
vector sent over the channel. The decoder tests the received vector y for containment in the 
code; if the test fails, it detects an error. Thus, the only case when the error is not detected 
occurs when the received vector y£ D \c. The probability to receive a vector y if the sent 
vector is c equals P(y|c) = (^r) ' (1 — p) n ~ dlst (y> c ). Then the probability of undetected 
error for the code D equals 



Pue(D, P ) = ^-S2 £ P(y|c) = -^£ £ ( ^_)-Hy,c) (1 _ prdl st( y ,c) 
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The theory of error detection for classical codes is surveyed in M . 

Observe that this (classical) concept consists of two parts, the definition of the error event 
and a way, (1), to compute its probability. In the quantum case both parts are not nearly 
as obvious. One can consider error detection with stabilizer codes; this produces a definition 
similar to the above one. It is not difficult to show that the probability of undetected error in 
this case can be computed via the weight enumerators of quantum codes in a way analogous 
to (1). This definition can be generalized in two different ways. First, it is desirable to extend 
it to cover all quantum codes; second and more importantly, from the physical point of view if 
the angle between the received vector and the transmitted vector is very small, it is natural to 
assume that no error has occurred since the measurement of the state of the system produces 
almost the same state as the transmitted one. On the other hand, if the received vector is a 
valid code point but is orthogonal to the transmitted one, then with probability one the error is 
not detected. Hence it is natural to define the undetected error event as the average probability 
that the state received from the channel is orthogonal to the transmitted one. 

A further generalization, which is physically perhaps the most important, involves trans- 
mission of completely entangled states and studying the probability of preserving the original 
entanglement. 



The goal of our paper in the first part is to prove that in all the situations mentioned, 
one can introduce a consistent definition of the undetected error event, and that the actual 
functional on the quantum code accounting for this event is essentially the same. This enables 
us in the second part || to prove that there exist quantum codes with exponentially falling 
probability of undetected error and to derive bounds on this exponent. 

The first part of the paper is organized as follows. In Section |2| we recall some notions of the 
theory of quantum codes, most importantly, the weight enumerators. In Section |3] we define and 
compute the probability of undetected error for stabilizer and unrestricted codes. The answers 
differ by a constant factor; hence their dependence on the code is the same in both cases. In 
Section ||| we study the event of undetected error in the case of completely entangled states 
and again arrive to the same functional as in the previous sections. Section |B| is devoted to the 
study of polynomial invariants of quantum codes. We show that "unitary weight enumerators" 
of pq] are known in classical coding theory under the name of binomial moments of the distance 



distribution. They were introduced by Mac Williams in [12]. They were studied extensively in 
[]l]] (see also related works referenced there), partly because they are convenient for bounding 
below the probability of undetected error. 

2 Quantum Error Correcting Codes 

In this section we review basic facts of the theory of quantum codes relevant to our study, 
focusing on the error process in the depolarizing channel and weight enumerators of codes. 

We begin by setting up basic linear-algebraic notation and reminding the reader of ele- 
mentary quantum-mechanical operations performed on state vectors. General sources for the 
relevant aspects of quantum theory are books |8[, |13|| ; a treatment that highlights the context 
of information transmission is given in |fL5||. Below by bold letters (v, w, . . . ) we denote complex 
column vectors. For a given vector v we denote by v* its conjugate transpose. Let TL n = C 2 
denote the complex 2 n -dimensional space. Let us fix an orthonormal basis in this space; denote 
it by vi , . . . , V2« . Observe that once the basis is fixed, a basis vector can be referred to by its 
number, i. This is employed when the states and operations are written in the Dirac notation, 
used in much of the physics literature on quantum codes. The correspondence is established as 
follows: 

Vi<r+\i), v* <-> (i\, v* Bvi <-> (i\B\i) , 

where B is a 2" x 2" matrix. This correspondence is discussed from the physcal perspecttive 
in []13|; mathematically oriented readers might enjoy the discussion in ]10| ]. 

Let V and W be two subspaces of TL n and A and B be linear operators on V and W, 
respectively. Consider the linear operator C = A®B on V® W. Let {vi, V2, . . . } be an 
orthonormal basis of V and {wi,W2, . . . } an orthonormal basis of W. The partial trace of C 
over V by definition equals 

Tr v (C) := Ti(A)B = £(v?Av 4 )B. 

i 

Here Tr(A) = Try (^4) = 2j( v I^ v i) is the trace of A. Similarly, one can define the partial trace 
over W as follows: 

Tr w (C) := ATr(B) = ^^w^), 



Obviously, 

Thw(C) = Tr v (Tr w (C)). (2) 

Any linear operator C on V (8> W can be written in the form 

j 

where Aj and Bj are operators on V and W, respectively. This representation is generally not 
unique. However, it is not difficult to check that the partial trace 

Trv(C) = £ T ^j)Bj, TV W = ]T Aj Tr(Bj). (3) 

3 3 

is a well-defined function. (In fact, it is possible to give an invariant definition of the trace [14, 



p.130].) 

Let us proceed to the definition of quantum codes and operations. A qubit is a two dimen- 
sional Hilbert space. A state (more precisely, a pure state) of a qubit is a unit vector v in the 
Hilbert space C 2 . Physically, a qubit occurs as a spin one-half particle, for example. Qubits are 
combined to form larger systems by taking the tensor products of the Hilbert spaces. The state 
of n qubits is therefore described by a unit vector in (C 2 )® n (nth tensor power). A quantum 
code Q is a linear subspace of 7i n . We denote a code of dimension K by Q((n, K)). 

Remark. There is a subtle point about this definition. Namely, although Q is a linear 
subspace, proportional vectors account for one and the same state, so in effect we only deal 
with vectors of unit norm. (Below we write u oc v to indicate that the vectors u and v are 
collinear.) Therefore, Q can be thought of as the projective space PC^ _1 . However, defining 
Q in this way poses problems for tensoring it with other subspaces (systems). Therefore, we 
prefer the definition given above. On the other hand, averaging over the code, we integrate only 
over the vectors in Q of unit norm. This convention, adopted with some abuse of notation, is 
valid for the entire paper. 

The number Rq = -^ — is called the rate of Q. A code vector (a state of n qubits) "sent" 
over a quantum channel is subjected to an error process that can alter the amplitude and/or 
the phase of some of the qubits. At the receiving end the decoder attempts to recover the state 
sent. The decoder is a quantum computing device; so it can perform only unitary rotations and 
measurements. Let us proceed to describing them in more detail. 

Let v e TC n , ||v|| = 1. A unitary rotation U of the state v 6 TC n produces a state Uv. A 
measurement of v is a probabilistic operation performed with respect to a set of projections 
Pi, P2, ■ ■ ■ ,Pt that satisfy the properties 

Pi + P 2 + . . . + Pt = id. 

Measuring v with respect to P\,P2,... ,Pt results in projecting it by one of the operators 
Pi, 1 < i < t. The resulting state has the form 

yv*PjV 



for some i. The probability for the ith projection is v*Pjv. Upon performing the measurement 
we obtain the state (||) and observe the number i. 

Suppose a code vector v € Q is sent through the depolarizing channel. The outcome of the 
transmission can be written as Ev, where the error operator E has the form 



E = Tl <g> T2 <S> ■ ■ ■ <8> T n , 

and each n is either I2 or one of the following (Pauli) matrices: 
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The weight of the error E is the number of nonidentity matrices in the expansion (|5|). Matrices 
of the form (|5|) either commute or anticommute: 

-Ei E? = uiEoE 



2-ca- 



The choice of the sign is determined by either geometric || or algebraic considerations ||] ; see 
( |13|) below. Note also that error operators are "trace-orthogonal": 

Ti[EiEj) = 2 n <% (7) 

Throughout the paper we assume that the channel is symmetric, i.e., the probability 

Prfo) 



\ -p, n = h 

n £ {o- x ,cr y ,a z }, 



p 

. 3' 



where < p < 3/4. Different qubits are subjected to the error process independently. Therefore 
we have 



Pr 



^wt(E) /n n _ wt{E) 



m = S) l '(i- P r 



(8) 



Similarly to classical codes one can introduce the weight enumerator B{y) = Y^=o ^iU l OI " a 
quantum ((n, if)) code [18]. By definition, 



B t :=^- 2 E Tr 2 (PP). 



if 2 



(9) 



wt(E)=i 



Here P is the orthogonal projection on Q. Another weight enumerator associated with Q ||T 
is given by B ± (y) = J2t=o B^y 1 , where 



B 



1 
K 



E Tr (EPEP). 



(10) 



wt(E)=i 



Some reflection shows that both Bi and B^ are real. Indeed, by (pi) error operators E are 
Hermitian. Of course, P is also Hermitian, and so is PEP. Eigenvalues of Hermitian operators 
are real; therefore, both Tr(EP) and Tr(EPEP) are real. 

The following two theorems account for the role of the weight enumerators in the theory of 
quantum codes. 



Theorem 1 p] 



B(y) = ^ w B(l + 3y,l-y) (11) 



Theorem 2 [18] Let Q be a quantum code with weight distributions B{ and B^ . Then 



i) B = Bfr = 1 and Bf > Bt > (1 < i < n) 

ii) the minimum distance of Q equals t + 1, where t is the largest integer such that Bi = 
B±,0<i <t. 

These theorems are readily seen to be analogous to the MacWilliams transform and sur- 
rounding results in classical coding theory. A particular case of classical codes is formed by 
additive (or group) codes |], p. 84]. An additive code C is a subgroup of the additive group 
of Z g . For such codes it is possible to define a dual code C , and the Hamming weight enu- 
merators of C and C 1 - are connected by the MacWilliams equation (|i~l|). A similar concept, 
stabilizer codes, was introduced in quantum coding theory in Q, Jj]]. 

To define stabilizer codes, let us provide more details on error operators (||). A natural 
framework to study them is that of orthogonal geometry and the theory of finite groups ||. 
Matrices of the form aE, where a 2 = ±1 and E is given by (||), form a group, say £ n , isomorphic 
to an extraspecial 2-group of order 2 2n+2 PJ. Its center is Z = {±/2 n ,±^2 n }- Let £ n = £ n /Z 
be the quotient group, formed by matrices of the form (||) with constant factors disregarded. 
In what follows an error operator E can be either from £ n or £ n . If E £ £ n then we will assume 
that E is of the form (|5|), that is we take E as a coset representative. 

Definition. A quantum code Q is called a stabilizer code if there exists an Abelian subgroup 
S of the group £ n such that Q is an eigenspace of S with eigenvalue 1. In other words, a code 
Q is a stabilizer code if 

Q = {veH n |V EeS £v = v}. 

If the order of S is 2 n ~ k , then dim(Q) = 2 k ||, @, and so the dimension of a stabilizer code 
always equals an integer power of 2. A stabilizer code of length n and dimension 2 is denoted 

by Q[[n,k]}. 

Let F := GF{4) = {0, 1,lo,lo 2 }. In the remaining part of this section we explain a way to 
associate to a quantum stabilizer code two additive codes C and C over F in such a manner 
that B(x, y) and B^(x, y) are their respective weight enumerators. This connection is the main 
discovery of ||]. It is also of key importance for both parts of our study. 

We begin by establishing a bijection (f> between the set of matrices {l2,o~ x ,a z ,a y } and the 
elements of F as follows: 

4>(I 2 ) = 0, <j)(a x ) = 1, 4>(a z ) = co, 4>(a y ) = to 2 , 

where a; is a primitive element of F. This bijection naturally extends to a bijection c/) n between 
£ n and F n and defines a mapping <p n : £ n — ► F n as a composition of (j) n and the factorization 



mapping. For a £ F let a be its conjugate under the action of the Galois group. In particular, 
this action transposes u> and Co = lo 2 . Let * be the inner product in F n defined by 

u * v = tr( Y^ Vi^i) = ^2( v i u l + v l u i)i ( 12 ) 

i i 

where tr(:r) := x + x 2 is the trace from F to GF{2). Key properties of error operators from £ n 
in this context are the following Q, ||]: 

El E 2 = l E2El e i* e ^=°' (13) 

[-£ 2 #i ei*e 2 /0 

where ej = cj) n (Ei), and 

E 1 E 2 = E 3 =^ ei + e 2 = e 3 . 

Thus, 4> n establishes a group isomorphism between £ n and the additive group of F n ; so £ n is 
elementary Abelian of order 2 2n , and 4> n is a homomorphism. Another useful property of error 
operators is 

Tr(E 1 E 2 E 1 E 2 ) = (-l) ei * e2 Tr(E l E 2 E 2 E l ) = (-1)^2™. (14) 

Now let S < £ n be a subgroup. Let 

C:={ceF"|c = <f) n (E) for some E e S}. 
Obviously, C is an additive subgroup in the group F n . Let 

C ± = {c£F n \ V c , ec (c * c') = 0.} 
Letting c range over the entire C 1 - defines a subgroup S 1 - = {^~ 1 (c)} < £ n . By definitions 



of C, C , and (|13| ) we see that any matrices E £ S, E' G 5 commute. Note also that S is 
Abelian; by (|i~3|) this implies that c' * c" = for any vectors c', c" G C; hence 

and therefore, 5 < S 1 - and 5 < S . 

The following theorem completes our task. 



Theorem 3 W% Let Q be a stabilizer code, B(x,y) and B (x,y) its weight enumerators, and 
let C and C be the codes over F n defined above. Then the Hamming weight enumerator of C 
(resp., C- 1 ) is B(x,y) (B ± (x,y)). 

The proof relies on the following lemma whose short proof we include for completeness. 



Lemma 4 [17] The orthogonal projection P on Q can be written in the following form: 

1 



On—k 



E E - ( 15 ) 



Proof A linear operator is a projection if and only if it is idempotent. A projection is orthogonal 
if and only it is self-adjoint. Self-adjointness of P follows from (|6|) and (|5[). P is also idempotent. 

? 2 = i-^) 2 E # E e = (^-) 2 E * n - k E = p- 

E'&s Ees EeS 

Let us compute the dimension of the space {Pv | v G TL}. From (J5|) we see that Tr(P) = 
unless E = 72™ • Hence Tr(P) = 2 k = dimQ. Finally, for any v G Q we have Ev = v by 
definition of Q; hence also Pv = v. □ 

The proof of Theorem (||) is now accomplished by invoking ( |l4| ) and (0) . 



3 Quantum undetected error 

In this section we study error detection with quantum codes used over the depolarizing channel. 
We examine a number of possible definitions of undetected error and establish their equivalence. 

For a quantum code Q we denote by Q 1 - its orthogonal code under the standard Hermitian 
inner product on 7i n . The orthogonal projection on Q is denoted by P . 

Let us first describe the quantum error detection protocol. The transmitted vector vGQis 
corrupted in the channel by the action of an error operator E G £ n . The received vector has the 
form w = Ev. In order to detect the error we perform the measurement of w with respect to 
the projections (P,P). If the result of this measurement is contained Q we detect an error. 
Otherwise we assume that there is no error and that the code vector Pw was transmitted. 

A. stabilizer codes Let us begin with the easier case of stabilizer codes. Suppose S is a 
stabilizer group of a code Q. Below by C we denote the code obtained by evaluating (f> n (E) on 

all E G S. 

Let E G £ n and v G Q. Suppose z is a result of the measurement of w = Ev with respect 
to the system (P,P). Let us introduce the indicator function 5(E, v) as follows: 

r/rn x f 1 z G Q and z gt v 
o(E,v) = < 

[0 otherwise. 

Note that if w G Q then measuring it with respect to the system (P, P ) with probability 1 we 
obtain a code vector. As shown below, for stabilizer codes the vector w is contained either in 
Q or in Q ; hence in this case 5 is deterministic. 

Definition. Let Q be a stabilizer code. Its probability of undetected error in a depolarizing 
channel with error probability p is defined as 

P ( uKQ,P) = E Pr (^) / 5(E,v)du(v), (16) 

where v(-) is the normalized uniform measure on Q. 

To compute the probability of undetected error let us first specify the error event, i.e., 
examine the case S(E,v) = 1 in more detail. Consider a possible effect of error operators 



on vectors from Q. For an error operator E let e = <p n {E). There are the following three 
possibilities. 

1) If E G S then Ev oc v for any v G Q; so this error has no effect on the transmission. 

2) Let E G" S . This means that e G" C 1 - and so e * e' ^ for some e' € C. Hence 
EE' = —E'E and E'Ev = —EE'v = —Ev. Thus Ev is contained in the eigenspace of E' with 
eigenvalue — 1. Since .E is a Hermitian operator, its eigenspaces are pairwise orthogonal. Thus 

w*Bv = 

for any w G Q. In other words w G Q^, and making measurement with respect to (P, P^) we 
will detect the presence of error with probability 1. 

3) Now let BeS^S, i.e., e e C ± \C. Then e * e' = for all e' G C. Hence by (0) 
EE' = E'E for all E' G <S, and so E'Ev = EE'v. Hence Ev is an eigenvector of E' with 
eigenvalue 1 and therefore Ev G Q. However, E ^ S; so Ev gt v for some v G Q. Thus, E acts 
on Q by rotation, possibly with some invariant directions. 

Note that in contrast with classical linear codes there are errors whose detection depends 
on the transmitted vector. 

The above information enables us to compute the probability of undetected error. 
Proposition 5 

Pi s J(Q,p)= E Pr (£) ( 17 ) 

egs- l \s 
= J2(B^-B i )^J(l-pr- i (18) 

Proof. The proof follows the standard proof of (1), relying upon @ and the definition of the 
weight distributions {Bi} and {B^~}. 

Suppose E G S 1 - \ S and e = 4> n (E). Let 

Q(E) -={v£Q\Ev = v}. 

Let us construct the code C e by adjoining e to the basis of C. Since e * e = for any e G C 2 , 
the new code is self-orthogonal. Hence we can construct an [[n, k — 1]] quantum stabilizer code 
Q{E) associated with C e . Obviously Q{E) C Q and any v G Q{E) is stabilized by E. It also 
follows from the definition of stabilizer codes that all vectors of C 2 stabilized by both S and 
E belong to Q(E). Hence Q(E) = Q(E). Therefore, the size of Q(E) does not depend on E, 
and for every E the set Q(E) is a K — 1-dimensional subspace. Therefore u(Q(E)) = 0. This 
implies our claim since 

P£HQ,P) = Y, F *( E ) f S(E,v)du(v) 

= E Pt(E)v(Q\Q(E))= E Pr(£MQ) 

Ees^s EeS->-\S 

= E Pr (^)- 

se5 x \5 



This proves (17). To prove the final part, note that \{E : wt(E) = i,E G S 1 - \ S}\ = B^ — Bi. 
□ 

This proposition shows that if, speaking loosely, we assume that the error E is not detectable 
whenever E G S \ S, ignoring vectors from Q, the probability of undetected error will be still 
given by (17)-(18). This might motivate another definition of undetected error for stabilizer 
codes, namely, assuming the uniform distribution on Q, let us call an error E undetectable if 
E G S 1 \ S. 

B. nonstabilizer codes. In this part we consider a more general situation of Q an arbitrary 
quantum code. Moreover, we take a viewpoint more natural in the physical perspective, namely 
if a "received" vector is close to the transmitted vector then we assume that little error has 
occurred despite the fact that formally they are not equal. This assumption is justified since a 
physical measurement in this case will not exhibit any difference between the transmitted and 
the received vectors. 

More specifically, let v G Q and w = Ev for some error operator E. If as a result of the 
measurement with respect to (P, P -1 ), w is projected on Q- 1 , the error is detected. Let us 
examine more closely the situation of w projected on Q. Let 

Pw PPv , x 

z = = - 19 

Vw*Pw \/v*EPEv 

be this projection, which occurs with probability v*EPEv. A natural proximity measure of 
z and v is the absolute value of the angle Z(z,v) between them. If z is sufficiently close to 
v, we assume that no error has occurred. By way of thought experiment suppose that for 
a given code vector v we measure z with respect to the system (vv*,P — vv*). Denote the 
result of this measurement by y z . In other words, we represent Q as a direct sum of v and a 
K — 1-dimensional orthogonal subspace Q v and create a pair of projections, on the line given 
by v and on Q v . The probability that after this measurement z G Q projects on v equals 
z*vv*z = cos 2 Z(z,v); the probability of the complementary event is z*(P — vv*)z. In the 
former case we assume that the error has no effect on transmission; in the latter that the error 
is undetectable. 

The overall probability that if v £ Q is subjected to an error operator E, then y z is a code 
vector orthogonal to v equals ||(I — vv*)PPv|| 2 . Indeed, it equals 

Pr{z G Q}Pr{y z _L v|z G Q} = (v*PPPv)(z*(P - vv*)z) 

^ v *PP(/ 2 n - vv*)PPv 

( = v*EP(I 2 n - vv*)(J 2 « - vv*)PPv 

= 11(7- vv*)PPv|| 2 , 



where (a) follows upon substituting z from (19) and recalling that P is a projection on Q and 



(b) uses the fact that I 2 " — vv* is a projection. Concluding, we arrive at the following general 
definition of the probability of undetected error. 



10 



Definition. Let Q be a quantum code used over a depolarizing channel with error probability 
p. Then 

P£\Q,P) =E Pr (£) / ll(/-w*)P^v|| 2 ^(v), (20) 

E J veQ 

where v is a normalized uniform measure on Q. 

This definition is more general than the one given for stabilizer codes. 

The main result that we prove regarding Pue {Q,p) is given in the following theorem, which 
shows that Pue (Q,p) differs from Put (Q,p) only by a constant factor. 

Theorem 6 Let Q be an ((n,K)) quantum code with weight distributions B{ and B^ , < i < 
n. Then 



^ ) (Q,p) = ZTT E(^ ± -^)(|) (i-p) n " J - (2i) 



i=0 

Proof. The proof is accomplished by computing the integral in (pO). It is clear that 

/ \\{I -vv*)PEv\\ 2 dv{v) 
JveQ 

= [ v*EP(I - vv*)(J - vv*)PEvdv(v) 
JveQ 

= f v*EPEvdv(v) - f (v* EPv)(v* PEv)dv(v). (22) 

JveQ ' JveQ 

To compute the integrals in (22) we need the following lemmas proved in the appendix. 

Lemma 7 Let Q be an (n,K) quantum code. Let P be the orthogonal projection on Q, and let 
//(v) be the normalized uniform measure on Q. Then 

f P 

/ vv*d^(v) = — 

JveQ & 

Lemma 8 

- ' (23) 



LeQ W * ® W *^ (V) = K(K + 1) (? Vj< ® ^ + V ^* ® Vj V * ) 



Using these lemmas, let us compute the integrals in (22). By Lemma we have for the first 
term 



/ v*EPEvdfi{v) = f Ti(EPEvv*)dfi(v) 
JveQ JveQ 

Tr (EPE f vv*d//(v) 
v JveQ 



^Tt(EPEP), (24) 



11 



where (a) is obtained by replacing a scalar by its trace. 

Let us compute the second term in (22). The code Q is a if -dimensional linear space; so it 
is isomorphic to the complex space C . The following calculations are simplified by performing 
them on C A instead of Q. 

Note that for any operators 

Ti(ABC) = Tr(BCA). (25) 

Therefore, we can rewrite our integral as 

/ \*EPw*PEvd/j.(\) = I Tr(EPw* PEvv*)dfj,(v). (26) 

Let L n be the space of linear operators on TC n . Consider the bilinear functional 

$i : L n x L n -> C 

(A,B) ^Tt(CAC*B). 

By the definition of the tensor product [J10| , there exists a (universal) bilinear map t : L n x L n — > 
L n <g> L n such that $i = $2 ° t, where <3?2 is a linear functional defined by 

$2 : L n (g> L n -> Tin 

A®B^Tt(CAC*B). 

Now we have 

/ v*EPw*PEvd(j,(v)=[ $ 2 (w*(g>vv*)d/u(v) 

Jvec K JveC K 

(b) 

^2 1 / vv 

V JvGC K 

K 



$2(1 vv* (g> w*d/j,(v) 

1 K 

^ ' »,f=l 



(c) 

~ K(K + 1) I ^ ^ Pv * v * PS E Vi v ■ ) + E TriEPvt^PEvjyt) 

~ K(K + l) [^ Tr{EPVlV * iPEP) + S Tr ( v *^ Pv *) Tr(v*P^ Vj )] 
(f) 



l—j [ Tr(PEPEP J2 Vi<) + Tr(EP £ v,v*) Tr(P£ £ v,-v*; 
(g) 1 



*(* + !) 



K{K + Tr(PEPE) + Tr(EP) 2 \ (27) 

Here (a) follows by ( |2"6| ) and the definition of $2) (b) holds true by linearity of <3?2i (c) follows 
by Lemma [|, (d) is by definition of §2, (e) uses the fact that P = J2i v^v*, (f) uses (|25|) and 
the additivity of trace, and (g) follows by (^5|) and the fact that P is idempotent. 

Substitution of (24) and (27) in (22) gives 

1^ \\(I-vv*)PEv\\ 2 dv(v) =J^- l [±Tr(EPEP) - ± Tv(EP) 2 ] ; 
by (|9|)-(^0|) this proves the theorem. □ 
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4 Composite systems 

In this section we study the most general problem setting for quantum error detection. The 
general idea is to take into account not only the error process but also the relationship of the 
current code vector with states of other quantum systems. More specifically, the qubits of the 
current vector can be entangled with other qubits that may not even take part in the trans- 
mission. We would like to study error detection that takes into account not only the error 
process but also this entanglement. The overall goal is to evaluate how well the original entan- 
glement is preserved under the action of errors. Somewhat surprisingly, though this definition 
of undetected error is more broad, the actual functional is again the same as studied above. 

This problem has no direct analogy in classical information transmission where typically 
one can study the effect of the error process on the transmitted vector without considering the 
influence on it of other parts of the message. 

Let Q be a quantum code that is a part of a combined system QR = Q (g) R, where R is 
another X-dimensional subspace. A generic element of QR can be written as J2i=i v « ® w i> 
where Vj and Wj are basis vectors of space Q and R respectively. Upon normalization we obtain 
a completely entangled state of the composite system, 



>0 * = vT?g v ' 



W; 



Let us assume that we transmit or store in quantum memory only qubits of Q. These qubits are 
subjected to the error process described above; so qubits of R remain error-free. The "received" 
state has the form 

1 K 

bQR = —j= ^{Eq ® I R )\i <g> Wj, 



where En is an error operator on Q and Ir is the identity operator on R. At the "receiving 
end" we again apply the same decoding, namely, measure the state with respect to (Pq,Pq). 
As above, the error is not detected if after this measurement we obtain a code vector that 
is orthogonal to the transmitted vector v. Therefore, in analogy with ( p0| ) let us define the 
probability of undetected error for the case of composite systems as follows 

i^ ) (Q,p):=E Pr (^)IK J «® J «- b O« b QH)( P Q® I B)( £; Q® J «) b Q«l| 2 ( 28 ) 

The main result of this section is given by the following theorem. 

Theorem 9 Suppose we transmit completely entangled states bg/j over a depolarizing channel 
with error probability p. Then 

where Bi , B^~ are the weight enumerators of the code Q . 
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Proof. Let us compute the summation term in (pq). We have 

\\(I R ® I R -b QR b* QR )(P Q ® I R )(E Q ® / /? )b QiJ || 2 

= b Qij[(^Q ® ^?)(^q ® ^)(^q ® Ir ~ bQ R h* QR )(I Q ®I R - bg fl b^) 

x (/q ® Jr - b Qfl b^)(P Q ® Ir)CEq ® Ik)] b Qi? 

= h* QR {E Q ® I R ){P Q ® I fl )(^Q ® / R )b Q R 

- (b^(£ Q Ir)(Pq I R )b QR )(b* QR (P Q ® I H )(£ g Jfl)b Qfl ). (29) 

The first term in this expression can be computed as follows 

Vqr(Eq ® I r )(Pq ® Ir){E q ® I R )h QR 

( = } Tr QR ((E Q ® J fl )(P Q 8) Jr)(£q ® /R)b QK b^) 

( = } Tr Q ( Tr fi ((£ Q I R )(P Q Jr)^ / R )b QR b^) 

= Tr Q ( Tr# (— J^ EqPqEqw^ ® WjW* 

(d) 1 

(e) 1_ 
if 
1 

where (a) is obtained upon replacing a scalar by its trace and using (pop, (b) follows by property 
@ of partial traces, in (c) we substitute the definition of \)q R , in (d) we compute the trace 
over R and use (pq), and in (e) we convolve over the dumb indices i,j. 

Let us compute the second term in (29). Proceeding as above, we obtain 
(b* QR (E Q ® I R )(P Q I R )b QR )(b* QR (P Q ® Ir)(£q I R )b QR ) 

= TT QR ((E Q ® /jj^Pq Ifl)bQ B b^) 

x T>Qij[(^Q ® Ifl)(^Q ® //?)b QiJ b^] 
= ^2 ^q(EqPqPq) Tr Q (P Q E Q P Q ) 
= ^Tr Q (P Q E Q f. (32) 



/f Tr <2 ( 5Z E Q P Q E Q w i w *j ® wfwiwjwj) (30) 



Tr Q (E Q P Q E Q P Q ) , (:{l) 



Substitution of (31) and (32) in (29) together with (|S|)-(|lC|) completes the proof. D 

This concludes our main task for the first part of the paper. We have proved that there 
exists a consistent definition of the probability of undetected error for quantum codes that can 
be given in the general case under natural physical assumptions and in the case of stabilizer 
codes analogously to the classical error detection. The functional of undetected error on Q in 
all the cases considered is the same, up to a constant factor. Therefore, as in classical coding 
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theory, we can study performance of quantum codes under error detection. The most important 
question in this context is to prove that the probability of undetected error for the best possible 
codes falls exponentially with code length n and to exhibit specific bounds on this exponent. 
Namely, let 

P ue (n,K,p) = min P ue (Q,p). 

Q£nn 
dim(Q) = K 

The second part of the paper is devoted to the study of this function. We answer the main 
question in positive by proving the existence of families of stabilizer codes with exponential 
decline of P ue (K,n,p), and establish upper bounds on this function for all quantum codes. 



5 Quantum weight enumerators 

In this section we focus on different forms of quantum weight enumerators. We begin with 
a short digression on classical enumerators. Let D C F™ be a linear code, supp(x) = {e G 
{1, 2, . . . ,n}|x e ^ 0} the support of a vector x (so wt(x) = |supp(x)|) and supp(P) = 
U xe x)supp(x) for a subset T> £ F™. Let (-£>i,0 < i < n) be the weight distribution of D. 
The following fact, proved by Mac Williams |l2| , underlies the combinatorial duality of coding 
theory: 




J2 \ V \ (0<w<n). (33) 



T>CD 

supp('D) | <w 



In particular, this immediately implies the Mac Williams identities [12). Denoting B w = J2T=o B% 
we obtain 



' n—i \ 
•n—w) 



Binomial moments of the weight distribution of codes were studied extensively in Q and some 
other related works (see [Q] for a discussion and complete bibliography). One of the reasons 
for this interest is that while any particular weight component Bi can be small relative to the 
code size, the numbers Bi can not. The probability of undetected error (1) can be written in 
the form 

pucp) = ±(b w - d(i - -^p^r 

10 = 1 " " 

thus lower bounds on B w are helpful for bounding P ue (C,p) below. 



For a quantum code Q one can generalize definition ( |33| ) by looking at error operators E 
whose supports are of restricted size, the support supp(-E) being the subset {e € {1, 2, ... , n}|r e ^ 
12}- Then we arrive at the following weight enumerators for Q: 

B w = ^ £ Tt^EP) 

ESSn 

I SUpp(B) I <TJJ 

Bi = ^ E Tr(EPEP). 

Eee„ 

I supp(.E) \<w 
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i A 

p«(0,p) = E(Bi-^)(f )( 1 -y) 



n—i 



The generating functions of these numbers, £>(x, y) and £> (x, y) were studied in [16| and called 
unitary weight enumerators. As in (^), it is immediate that 

B w = ± Bl ( n - l \ Bi = ±Bt( n - i \ 
to \ n - w ) to V n "W 

which is a result in [16]. The Mac Williams equation ( |l"l| ) also follows immediately by the 
original proof in [|E^]. Also, 

B(x, y) = B(x - y, y) B ± (x, y) = B ± (x- y, y); 

thus the probability of undetected error equals 

i=0 

So to bound P U e(Q,p) below we could first derive lower estimates on (Bi — B^) following the 
ideas of m. However, in part 2 we choose to work with the functions (18), ( |2l| ) as a whole. 

A Appendix 

We precede the proofs of Lemmas p] and M with two other useful facts. In the proofs below we 
repeatedly interchange the order of integration. Obviously, all the necessary assumptions on 
the measures (cr-additivity, completeness) for the Fubini theorem to hold true are in place. 

Lemma 10 Let G be a compact group and it be a unitary representation with operators acting 
on a linear space W . Let fi be the Haar measure on G. Then 



P ■= K g dn{g) 

JaeG 



IgeG 
is an orthogonal projection on W . 

PROOF. It suffices to show that P 2 = P and P* = P. 

52 



P = K h Tr g dfi(g)dfi(h) 

JheG Jg&G 



^hgdn(g)dfi(h) = / TT g dfi(g)dfi(h) 

heGJgeG JheGJgeG 

Since n is a unitary representation, ir* = vr„-i. Hence 

P* = ir*dfj,(g) = / Tc g -idfj,(g) = / n g dfj,(g) = P. 

JgeG JgeG J g£G 



□ 



As in the proof of Theorem || in the following lemmas we perform calculations in C K instead 
of Q. Define the inner product of matrices A and B as follows 

(A,B) =Tr (A* B). (34) 
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Lemma 11 Let L be the space of linear operators on C K and U be a unitary matrix. Then the 
equality 

A ® B = U* AU (g> U* BU (A,BeL) 

holds true for any unitary matrix U if and only if A®B is contained in the subspace W C L®L 
generated by 

K K 

I := I (g) I = V] v i v i ® v i v j an d J = z2 v * v j ® v i v i • 

Proof. The proof is at times sketchy, however we only omit routine calculations. Let it be 
the representation of U(K) acting on L L as follows 

nu(A ® B) = U*AU ® U*BU. 

We need to prove that iru(A (g> 5) = ,4 5 for all 17 € U(K) if and only if A (g) B E W. 

Let us begin with the "if" part. It suffices to prove that ttjj acts identically on I and J. The 
first of these is obvious. For the second, let us introduce the canonical isomorphism 

(f) : L®L^L®L 

A®B \-> B (g> A. 

Let us compute 7rj/(J) as follows: 

nu (J) = J2 U*ViV*U U*vjv*U 

y 

* i 

j 
53 Viv| (8)Vj v*. 

Let us prove the "only if" part. Consider the group G = DSk where D is the group of 
all diagonal matrices with diagonal elements from the set {±l,±i} and Sk is the symmetric 
group. Clearly, G < U(K). Let 

|G| gee 
By Lemma [H], \ is an orthogonal projection on G. We will prove that the dimension of its 
image is 3. After that we will present an operator w 6 imx that is invariant under the action 
of G but is not fixed by U(K). This will imply that the dimension of the subspace of L <S> L 
fixed by U(K) equals 2; hence by the above this subspace is W. 
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Let us find dimim(x) = Tr(x) = u*xuj, where (iij, i = 1, . . . , K A ) is any ortho-normal basis 
of L <g) L. For instance, take 

UjjfcJ = VjV* (g) VfcVj*. 

Then 

Tr ( igf E %) = ]gf E ( v * v ) ® v * v * > E %( v * v ! ® v * v f ))> 



■|G| 



\G\ / 



3GG ' ' i,j,k,l geG 

where (-, •) is the Hermitian inner product on L (g) L, 

{A®B,C®D) = Tr(AC* (g> BD*). 



(35) 



Consider an element g G G, g = diag(s)<r, where s G {±l,±i}^ and <r G Sk- It is easy to 
see that 

7T ff (VjV* VfcV*) = SjS*SjfcS*V (T (j)V* (j) V^V*^. 

Since the basis is orthogonal, we have 



*'-m£" 



i 



X] ( V i V j ® v fc v *> E ^(ViVj ® v fc v f) 



i,j,k,l 



geG(i,j,k,l) 



where G(i,j,k,l) < G is a subgroup formed by the elements g = diag(s) <r, where <r has fixed 
points i,j, k, I as fixed points. The inner product under the sum is nonzero only in the following 
three cases: 

a) i = j,k = l,i / k 

b) i = l,j = k,i / k 

c) i = j = k = I. 

It is easy to check that in each of these cases the sum equals \G\. Thus we have 



1 

\G\ 



Mtt^E^)^ 

geG 



To complete the proof notice that the element 

w = E v * v ** ® v * v ** G w 



is invariant under the action of G but not of U(K). Indeed 



K 



^ffOO = E V <?W V ki) ® V "(0 V *(i) = K 



i=l 



It is easy to check that V is not fixed under the action of the unitary K x K matrix 



U 



V2 



1 


1 


.. 


. " 


1 


-1 


.. 


. 








V2 .. 


. 








.. 


. v/2. 
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For instance for K = 2 we have 

r 1 1 



V2 



1 -1 



and it is straightforward to see that 



Uv^lU &> C/vivjC/ + Uv 2 v* 2 U ® Uv 2 v* 2 U 

= (vi + v 2 )(vi + v 2 )* <8) (vi + v 2 )(vi + v 2 )* 

+ (vi - v 2 )(vi - V 2 )* (g) (vi - v 2 )(vi - v 2 )* 
/ ViV^ (g) V]V* + v 2 v 2 v 2 v 2 . p. 

Proof of Lemma |7|. 

Let Wj and v, be orthonormal bases of Q and C , respectively. Let <fi : Q —> C be the 
natural isomorphism given by </>( w i) = Vj. Denote by p(U) the Haar measure on the unitary 
group U(K). Applying (j>, we can rewrite the integral in question as follows 

/ ww*rf/i(v) — > / vv*<i/i(v) 

JveQ Jvec K 

= f U*vv*Udfi(-v) 
Jvec K 

= f dp{U) f U*w*Ud[j,(y) 

Ju<=u(K) Jvec K 

= f U*vv*Udp(U), 

JUeU(K) 

where in the last expression v is an arbitrary fixed basis vector. We need to compute the last 
integral. 

Let it be a unitary representation of U(K) acting on the vector space of complex K x K 
matrices as follows 

7Ttf(w*) = C/*vv*C7. 

Consider the operator 

/ U*vv*Udp(U) = [ TCu(vv*)dp(U) 

JU£U(K) JU<=U(K) 



by Lemma 10 it is an orthogonal projection. Next we show that it projects on the one- 



dimensional subspace generated by Ik', then the last integral becomes easy to compute. 
Note that it is unitary with respect to the inner product (|34|). Indeed, 

(tt u (A),tt u (B)) = Tt(U*AUU*BU) = Tr(A*B) = (A,B). 

From the standard fact that any matrix A' can be represented in the form A' = UAW, where 
U and W are unitary matrices and A is diagonal, it is easy to see that the identity U*AU = A 
can hold for any unitary matrix U if and only if A = Ik (this is a particular case of Lemma 



11). Therefore 

U*vv*Udp(U) 

ueU(K) 
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equals the projection of v*v on the one-dimensional subspace of W generated by Ik- Recalling 
that the projection of a vector x on a vector y equals )*'y( y, we have 

JueU(K) {Ik, Ik) K K f^i 

Finally, since 

0- 1 (f;w < wj)=f;v 4 v? = p ) 

we are done. □ 

Proof of Lemma ||. 

Let L be the space of linear operators on C^. Similarly to the proof of Lemma it can be 
seen that 

/ vv*®vv>(v)=/ U*vv*U®U*vv*Udp{U), (36) 

Jv£C K JueU(K) 

where v is an arbitrary fixed basis vector. By Lemma [ll| the integral in (|36|) is a projector, 
with respect to the inner product ( |35| ) , on the subspace W C L (g) L generated by I and J. The 
lemma will be proved if we evaluate the projection on this subspace of u := vv* <g> vv*. 

To do this, we need an orthogonal basis of W (note that I and J are not orthogonal; indeed, 
it is easy to check that 

(T, I) = K 2 , (J, J) = K 2 , and (J, J) = K.) 

Therefore, let us consider the basis /, J, where 

Again it is easy to see that 

(J, J) = K 2 - 1 and (J, J) = 0. 

Thus, / and J form an orthogonal basis of W. To compute the projection of u on W, let us 
first project it on the basis directions: 



(J,u) = Tr(^(v i v*®v i v*)(vv*^vv > 
= ^Tr(v,v*vv*)Tr(v,v*vv*) 



h3 

= l; 

Thus, the projection of u on W" equals 

VV* ® W*dA*(v) = ^^/+ ^i^J = — ^ r(J+ J). 

This completes the proof. □ 
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